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: FACILITY: MATH LIBRARY 

¥ ABSTRACT: 

MTHSHSIN and MTHSHCOS are functions which return the Pde) ia point 
; sine or cosine value of their single precision floating point argu- 
3 ment (radians). The call is standard call-b arenes ence. 

3; MTHSHSIN RS got MTHSHCOS_R5 are special routines which are the same 


; as MTHSHSIN and MTHSHCOS” except a faster non-standard JSB call is 
; ; used with the argument in RO and no registers are saved. 


: : MTHSHSINCOS is a routine which returns the floating point sine and 
; cosine value of its single precision Hoos ioe 2 aot radian argument. 
; The call is standard call-by-reference. MTH S_R5 is a special 
; routine which is the same as MTHS$H HSINCOS. remy a faster non- 
; t are'seved. call is used with the argument in RO and no registers 
; are saved. 


 MTHSHSIND and MTH$HCOSD are fynct ions which return the Seating point 
sine or Bt ss value of their single © pres tenes floating point argu- 
: cont (de The call is standar ce l-by-reference. 
; MTHSHSIN H ae MT Wen OSD R5 are special routines which are the same 
; as MTHSHSIND and MTHSHCOSD” except a faster non-standard JSB call is 
; ; used with the argument in RO and no registers are saved. 


; : MTHSHSINCOSD is aro tine which returns the floating point sine and 
; cosine value of its single precision floating peret one's argument. 
; The call is standare call-by-reference. MTHSHSINCOSD_R5 is a Special 
; routine which is the same as MTHSHSINCOSD, except a faster non- 
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Edit History for Version 1 of MTH$HSINCOS 


1-001 = Adapted from MTHSGSINCOS version 1-001. JAW 21-Aug-1979 
1-002 - Improve argument reduction scheme for greater geguracy- Also 
change MOVH to MOVO in MTHSHCOS. JAW 28-Jan-1980 
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-001 = Changed W* modifiers to G* on external references RNH 09-Sep-81 
= Changed MTH$AL_4_OV_PI to MTHSAL_4_OV_PI_V RNH 29-Sep-81 

- Included check for A2=0 before decrementing exponent. RNH 02-Oct-81 

° anc ‘ges z gpeerrs on references to MTHSAL_4_OV_PI_V. 

-Uct= 
2-005 - Modified logic for converting reduced argument from tnteger to 
floating point format to avoid modifying exponent of a floating 

oint zero. RNH 21-0ct- 
2-006 - pee are | ct logic to eliminate bug reported in QAR 896. 
-Jan- 
2-007 - In after porere ine from a JSB call to CVT_TO_H, the code in CONVERT: 
did a SUBW #*X1D,16(R5), to change h_hi*2"29°to h_hi. However, this 
only valid if the contents of 16(R5) is not zero. The code has been 
patched to skip the above instruction if h_hi is zero. JCW 22-SEP-83 
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MTHSHSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:39: AX/VMS Macro V04-00 Pa 
2-007 ThSns in 8s = Radian arguments aes 7 94:38:50 ERTHRTL. SREIATHHSINCO.MAR; 1 - 
~SBTTL MTHSHSINCOS - Radian arguments 


FUNCTIONAL DESCRIPTION: 


The HSIN, HCOS and HSINCOS routines are based on octant reduction. Given an 
argument, x, it is written in the form 


x = [1e(2*pi) + I*(pi/4) + Y1, 
where 11 and I are intggers 0 =< I < 8 and 0 =< Y1 < pi/4. Since HSIN and 
od of it f 


WO OOO O00 000000090000 000008 SII 
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C 
9 
9 
9 
9 ; 
9 ; 
9 ; 
9 3 
9 ; 
9 3 
9 ; 
95 ; 
2 ; 
, ; HCOS have a per pi ollows that 
98 5 HSIN(x) = HSIN(I*(pi/4) + Y1)) and 

035 3 HCOS(x) = HCOS(I*(pi/4) + Y1)). 

95 ; Using the trigonometric identities for the sum and difference of two angles, 

Baap ; the following table can be generated: 

0950 96 ; If l= then HSIN(x) = and HCOS(x) = 

0950 97: jj eeeeem j= ewmmeemeseres = = commacceecece 

0950 98 ; 0 HSIN(Y1) co 

0950 99 ; 1 HCOS(pi/4-Y1) HSIN(pi/4-Y1) 

0950 400; ¢ HCOS “HSIN 

0950 401 ; HSIN(pi/4-Y1) “HCOS(pi/4-Y1) 

0950 402 ; 4 “HSIN(Y1) “HCOS(Y1) 

0950 403; 5 saereige mt 2 ata 

0950 404 ; 6 “-HCOS(Y1) HSIN(Y1) 

bae6 ret 3 7 “HSIN(pi/4-Y1) HCOS(pi/4-¥1) 

0950 407 ; Let Y be defined as Y = Y1 if I is even and Y = pi/4 - Y1, if I is odd, then 
950 408 ; each entry of the above table is of the for +/-HSIN(Y) or +/-HCOS(Y). Based 
950 409 ; on the above remarks, the HSIN, HCOS and HSINCOS routines process the input 
950 410 ; argument x, to obtain I and Y, and based on | selects a suitable polynomial 

0950 411; approximation, p(Y), to evaluate the desired fuction. 
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; Floa Point Sine, Cosine od Sincos 16-SEP-1984 AX/VMS o V04-00 Page 11 
rhs as = Radian arguments met ti 33: 4h YRTMRTL. 8 L.S SATHHSINCO. MAR; 1 ° (8) 
: INPUT PARAMETERS: 
: LONG =4 
x = 1*LONG ; x is input angle in radians 
sine = Scions : sine is HSIN(x) 
cosine = 3*LONG 3; cosine is HCOS(x) 


; Return sine and cosine of argument 


ENTRY MTHSHSINCOS, “M<R2, R3, R4, RS, R6,R7> 
MTHSFLAG_JACKET 
MOVAB G*MTHSSJACKET_HND, (FP) 


MOOD NAUE WN O OONO NEW 


Foy oy at ot ot ot at at at a 


6D O0000000'GF 9E 
3; set handler address to jacket 


= 


DWDEMOOODOONNMMNNVOQOOCOCOCOOOOOCOOCOCOCOOO O37 


Deal * ah val val ab ab Ab Ab Ah Ab Ab Ab db db dh db db db Ah Ab Ab db dh Ab hb dh dh hh A” teed 


3; handler 
138 
50 04 BC 70FD 4 MOVH 
000009F6"EF 16 434 JSB at AsnsiNcos R7 
08 BC 0 7DFD 435 MOVO RO, asine(AP) 
OC BC 54 7DFD 4 $ MOVO R4, @cosine (AP) 
04 6— 4 RET 
96F 438 
96F 439; 
96F 440 ; INPUT PARAMETERS: 
96F 441; 
00000004 Hs 2 LONG =4 
000008 O96F 444 arg = 2*LONG ; x is input angle in radians 
0000004 O96F 445 answer = 1*LONG 


96F 44 ~SBTTL MTHSHSIN 


; Return sine of argument 


003¢ -ENTRY MTHSHSIN, “M<R2, R3, R4, RS5> 
MTHSFLAG_JACKET 


MOVAB G*MTHSSJACKET_HND, (FP) 


71 


6D Q00000000'GF 9 
3 set handler address to jacket 
3; handler 


SN 


© oapcD000000 — — 
ann 
Wu 
Con 


50 8 BC 70FD 


MOVH 5h ht 
00000B3A'EF 16 


JSB MTHSHSIN_ ae’ 
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04 BC 50 7DFD 
04 


003¢ 


6D 00000000'GF 9€ 


50. 08 BC 70FD 
000008D ‘ff 16 
04 BC 0 7DFD 


04 


-SEP-19 1: MTHRTL.SR RCM MTHHSINCO.MAR; 1 
460 MOVO RO, @answer (AP) 
4 RET 
4 
4 
re SBTTL MTHSHCOS 
46 < 
a8e Return cosine of argument 
4 
a7 
t76 sENTRY MTHSHCOS, “M<R2, R3, R4, R5> 
474 MTHSFLAG_JACKET 


MOVAB G*MTHSSJACKET_HND, (FP) 
; set handler address to jacket 
3; handler 


5 

8 MOVH Borg cana e RO 
JSB MTHSHCOS_R5 

8 MOVO RO, @answer (AP) 

; RET 


(6 2 


=} 
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argument Xs to obtain I and Y, and based on | selects a suitable polynomial 
rs) 


approximation, p(Y), to evaluate the desired fuction. 


MTHRTL.SR STH HHSINCO.MAR; 1 ( 
9A 4 § -SBTTL MTHSHSINCOSD - Degrees 
9A 4 
Ae 
9A 486 t FUNCTIONAL DESCRIPTION: 
9A 4 : The HSIND, HCOSD and HSINCOSD routines are based on octant reduction. Given 
aA 2 3 an argument, x, it is written in the form 
gh a3 : x = 11#360 + 1#45 + Y1, 
9A re t where 11 and I are ingegers 0 =< I < 8 and 0 =< Y1 < 45. Since HSIND and 
9A 49% ; HCOSD have a period of 360 it follows that 
9A 496 : HSIND(x) = HSIND(I®45 + 1) and 
aA hth 3 HCOSD(x) = HCOSD(1*45 + Y1). 
O3A 499 ; Using the trigonometric identities for the sum and difference of two angles, 
Bea 289 ; the following table can be generated: 
O9A 308 ; If l= then HSIND(x) = and HCOSD(x) = 
O9A 04 : 0 HSIND(Y1) HCOSD(Y1) 
O9A 82 3 1 Heosptaseyt) HSIND(45-Y1) 
O9A 6; ¢ HCOSD(Y1) ee 
O9A 444 ; SIND (45-¥1) D(45-¥1) 
O9A 08 ; 4 “HSIND(Y1) rHeOSD<yt) 
O9A 509 ; 5 ate O50 eet) rHSIND (a5 Y1) 
O9A 219 ; 6 -HCOSD(Y1) HSIND(Y1) 
SoA gi) 3 7 “HS IND (45-Y1) Heos(45-¥1) 
09A3 218 ; Let Y be defined as Y = Y1 if I is even and Y = 45 - Y1, if I is odd, then 
O9A a 3 each entry of the above table is of the for +/-HSIN(Y) or +/=HCOS(Y). Based 
O9A 1? 3 on the above remarks, the HSIND, ACOSD and HSINCOSD routines process the input 
21$ : 
eis 
19 


‘ 
ae 
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SINCOSD = Degrees “SEP-1 1: 
80000008 A 1 sind = gol Ons 
C aA § cosd = 5*LONG 
OOFC RA ¢ -ENTRY MTHSHSINCOSD “M<R2, R3, R4, RS, R6, R7> 
9A 6 MTHSFLAG_JACKET 
6D Q0000000'GF 9E O9A MOVAB G*MTHSSJACKET_HND, (FP) 
9AC 3 set handler address to jacket 
9AC 3; handler 
9AC 
9AC 7 
of BC 70FD OQ9AC 8 MOVH ax¢ 
00000¢ ver 16 0981 9 JSB RTWBHS{NCOsD R7 
BC 0 7DFD 0987 0 MOVO RO, nd (APY 
OC BC 54 7DFD ger 1 MOVO R4, SE OC dCAP) 
04 a3e1 § RET 
9C 4 
09C 5 
09C 3$ 
003C¢ 395 23f -ENTRY MTHSHSIND “M<R2, R3, R4, RS> 
asc 335 MTHSFLAG_JACKET 
6D 00000000 ' GF 9E 09C4 MOVAB G*MTHSSJACKET_HND, (FP) 
09C¢B : set handler address to jacket 
09C¢B : handler 
09CB 
pecs 49 
50 08 BC 70FD 09CB 41 MOVH Borgtary, & 
00000CEB EF 16 0900 4g JSB HSIND_ RS 
04 BC 0 7DFD 09D6 4 MOVO RO, @answer (AP) 
3444 44 
04 0908 $2 RET 
09DC 46 
O9DC 47 
9DC 48 
003C Bape 23 -ENTRY MTHSHCOSD “M<R2, R3, R4, RS> 
O90E 51 MTHSFLAG_JACKET 
6D 00000000 ' GF 9E 9D MOVAB G*MTHSSJACKET_HND, (FP) 
9E 3 set handler address to jacket 
9E : handler 
9E 
9E 38 
98.8 C 70FD O9E 5 MOVH aarg (AP) RO 
060090 C fF 16 O9EA 34 JSB MTH Sucosb. R 
04 0 7DFD er 2? MOovO RO, @answer (AP) 
04 oF 3$ RET 
9F6 58 
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af eT -SBTTL MTHSHSINCOS_R7 
oF gi This routine computes the HSIN and HCOS of the G-format value of RO/R1. The 
OF 65 ; ponputet on is performed one of three ways depending on the size of the 
are rt: ; input argument, X: 
or 66 : ; 1) if ae < Ray then X is used directly in polynomial approximation 
9F6 68 ; 2) If pi/4 =< ixi < “Qepis4, then the subroutine REDUCE_MEDIUM is called 
Bare 4 3 to reduce the ergunent to an equivalent oroyment in radians, ve + ~ 
9F6 0; the octant, I, conta . the ior ougent. is then evaluated in 
4 1 3 po tynom als chosen ps a Tunctio to” compute HSIN(X) and ACOs <x) « 
OF (2 3 3) =< iX!, then the Scarourten’ REDUCE_LARGE is called to 
bar 7 : reduce “ye argument to an equivalent argument in cycles, Y¥ A the 
9F6 74 ; octant, contanining the argument. is then evaluated 
848 4 g 3 a4 pat OR chosen as a function of Sy to compute HSIN(X) one HCOS (xX). 
bore 4 MTHSHSINCOS_R7:: 
50 73FD O9F6 78 TSTA RO ; Check sign of x 
10 18 O9F9 579 BGEQ POS_SINCOS : 
ie oh 16 O9FB 580 JSB SINTOS 3 RO/R1 = HSIN(CIX!), R2/R3 = HCOS(X) 
0 72FD OA01 3! MNEGH RO, RO ; RO/R1 = HSIN(X) 
05 OA0S 5 ¢ RSB 
OA06 «58 
QA06 584 SINCOS: 
50 8000 8F AA QA06 585 IcW #*x8000, RO 3; RO/R1 = ix! 
AOB 586 POS_SINCOS: 
50 FSFO CF 71FD QAOB 587 CMPH H PI _OV_4, RO 3; Compare pi/4 with {Xi 
14 QA11 ee BGTR SMALC_SINCOS 3 No argument reduction. is necessary 
50 FSF8 CF 71FD OA13 89 CMPH H_9_PT_OV_4, RO ; Compare 9*pi/4 with ix! 
03 18 OAI9 590 BGEQ 1$ 
OODE 31 male 235 BRW LARGE_SINCOS ; Use special logic for !Xi > 9xpi/4 
pale 38 H 
pale Be : pis =< IK! < 9epis4 
SE 30 C2 OAIE 596 is: SUBL #48, SP ; Allocate 12 Longword on the stack 
OA21 597 ; _ for the reduced argument 
OOOOODAS'EF 16 oA 1 295 JSB REDUCE _MEDIUM ; Medium argument reduction routine 
Ae? 99 3; Y = reduced argument on the stack, 
Ae? 600 3 a = octant, R5 points to YHI. and 
Ae? 601 ; > Porase to YLO 
M3 64 7DFD OA27 O06 MOVO (R4), =(SP) 3 eng 
¢ 7DFD OA2B 4 MOVO (R5), =(SP) : Save Y 
02 Ack 604 PUSHR #*M<R2, R4, R5> > Save _R2, R4 
oppocers” 3 6 QA31 605 JSB M_CO ; RO/R3 = HCOS( 
0 7DFD OA37 6 § MOVO RO, 76(SP) 3 Save HCOS(X) on stack 
4 BA OASC 6 POPR #*M<R2, R4, RS> 3; Restore R2, R4, 
E 7DFD OA 608 MOVO (SP)+, (R5$ ; Restore YHI 
Sa E 7DFD OA4 60 MOvVO (SP)+, (R4) ; Restore YLO 
ogo gt "EF 16 QAGG 610 JSB MSN : RO/RS = HSIN(X) 
$5 7D0FD QA4C 61 MOVO 32(SP), R4 3 R4&/R7_ = HCOS(X) 
c? A51 $16 ADDL #48, SP : Pop 12 longwords off the stack 
0 nee 91? RSB 
Age gi ; Logic for small arguments. !X! < pi/4. 
ASS 616; 


SINCOS 


F729 


F8CA 


F7C5 


F896 


7E 8 7DFD 
4000 8F : H 
50 3FC7 BF BI 
50 E : SF 
E 56 for 
CF C 50 75FD 
4 3 7DFD 
E ERED 
CF ) 4 75FD 
54 8E 7DFD 
50 6E 64FD 
50 8E 60FD 
05 
os 6f0hllU 
54 08 70FD 
05 
54 6€ 6€ 65FD 
7E 54 7DFD 
CF OD 54 75FD 
54 6 7DFD 
6E 50 7DFD 
CF OD 54 75FD 
54 8E 7DFD 
50 ©66E 64FD 
50 6€ 60FD 
7E__«-6E 7DFD 
FFFFOIFF 8F CA 
1C AE D4 
6E 10 AE 63FD 
10 AE 20 AE 60FD 
6— 8€ ore 

02 1 
6— 87 
54 8E€ 62FD 
6—E 6€ 64FD 
6— 87 
6E 8 osFD 
54 E 62FD 
05 


—. 
23 
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DEAE AAAAAA AAA AAA AAA AAA AA AA AAA AA AAA OA AAA AAAAAASAAAAIAOSAASOSSAAS Wo 
SSD DDD DDD DDD. TTT BS ES BS BE BE EAN INIPPOPONPPO PUN SS So“ 


WO OONOUES WO 0Oon 
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SMALL_SINCOS: 
MOVO 


1$: ADDL 


2$: MULH3 
M 


3$: SUBH 
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t s 
Save space for HCOS(X) 
Save reduce 


VMS Macro - 
6-5 MTHRTL.SRCJMTHHSINCO.MAR; 1 (11) 2-( 

RO, -(SP) 3 (SP) = {Xi 
RO, #*x4000 3; Compare 1/2 with iX! 
es ; Sufficent overghang not available 
#*X3FC7, RO 3; Compare with 2*- 
1$ 3; No ee‘prans evaluation is needed 
(SP), (SP), RO 3; RO/ = X*X 
RO, -(S ; Put X*X on stack 
RO, #COSLENR1-1, COSTBR1; RO/R3 = HCOS(X) 
(SP), R4 : RG/R7 = Xex 
RO, (SP) 3 Save_HCOS(X) on stack 
R4, MSINLENR-1, SINTBR ; RO/R3 = pb a 
(SP)+, RG t R4/R7 = HCOS(X) 
(SP), RO 3 RO/R3 = X*qQ(X*2) 
(SP)+, RO ; RO/R35 = HSIN(X) 
#16, SP Clear stack 
#1.6, R4 RO/R3 = X, R4/R7 = 1.0 = HCOS(X) 
(SP), (SP), R4 3 R4/R7 = X*2 
R4, =-(SP 3; Save x*2 
R4, MCOSLENR2-1, COSTBR2; RO/R3 = Q(x*2) 
(SP), R4 t RG/R7 = X42 
RO, (SP) 3 Save _Q(X*2) on stack 
R4, MSINLENR-1, SINTBR ; RO/R3 = Q(X*2) 
(SP)+, RG : RG/R7 = P(X*2) 
(SP), RO 3 RO/R3 = X*P(X*2) 
(SP), RO ; RO/R3S = HSIN(X) 
(SP), -(SP) ; Save another copy of :X: 
#°XFFFFOIFF, 24(SP) : 
28(SP) : 16(SP) = XHI 
16(SP) ,(SP), ; (SP) = XLO 
32(SP), 16(SP) : 16(SP) = X + X 
(SP)+, (SP) ; (SP) = XLO®(X + XHI) = A2 
3$ 3; Check for A2 = 
(SP) : (SP) = A2/2 
(SP)+, R4 : R4/R7 = Q(Y¥*%2) = A2/2 
(SP), (SP) 3 (SP) = tl 
(SP) : (SP) = XHI*2/2 
#1, (SP) 3; (SP) = XHI*2/2 = 1 
(SP)+, R4 + R4/R7 = HCOS(X) 
#32, 3? ; Allocate 16 longwords on the stack 
P, R 3 ointed to by R5 
REDUCE_LARGE ; (R5) = YLO 

; (R4) = YHI Cin cycles) 

3; R2 = octant bit 


argument on stack 


C 8 
MTHSHSINCOS ; Floating Point Sine, Cosine and Sincos ys nb 01:39 
2-007 ATHSHSINCOS._R7 SEpaloee $138 
34 BB 0B14 ore $ 
opposer" F 16 QB16 75 L_cos 3 RO/R1 
0 7DFD 081 4 MOVO RO, 44(SP) 3 (SP) = “ue sath 
4 BA OB 67 POPR #*M<R2, RG, RSD :; Restore R2, R4 and R5 
10 A E 7DFD 678 MOVO (SP)+, 16(R5 ; Reduced argument on stack 
6 — 7DFD 67 MOovVO (SP)+, (R5) 3 £9 by R4 and R5 
eins > alll 16 § 680 JSB L_SIN 3 RO/R3™ : ASIN(X) 
4 f 7D0FD OB 681 MOVO a RG 3: R4/R7 = HCOS(X) 
ps4 4 3 6 § ADDL #32, SP ; Pop 8 lLongwords off the stack 
B 6 RSB 


:11 VAX/VMS Macro V04-00 Page 17 au 
6-SEP=19 :20 CMTHRTL.SRCIMTHHSINCO.MAR; 1 (11) 
PUSHR #*M<R2, R4, RSD ; Save Re R2, we “itt RS 
JSB : 
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MTHSHSINCOS Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:39:11 VAX/VMS Macro V04-00 Page 18 
2-007 ATHSHSIN. RS g- aati 94:38:30 EMTHRTL. SRCIMTHHSINCO.MARS 1 . (12) 
: ; 685 -SBTTL MTHSHSIN_RS 
B3A 3 : ; This routine computes the HSIN of the G-format value of RO/R1. The 
BSA 688 ; computation is performed one of three ways depending on the size of the 
B : 283 ; input argument, X 
B rs $3) : 1) if Xi < pi/4, then X is used directly in a polynomial approximation 
§ ) ° 
BSA $38 ; 2) If pi/& =< ixt < 9epi/4, then the subroutine REDUCE_MEDIUM is called 
OBSA 694 ; to reduce the argument to an equivalent argument in radians, Y, and 
OB3SA 695 ; the octant, I, containing the argument. Y is then evaluated in a 
68 A 036 $ pol ynoayal chosen as a function of I to compute HSIN(X). 
BSA 697 ; 3) If 9epi/4 =< IX!, then the subroutine REDUCE_LARGE is called to 
OB3SA 698 ; reduce the argument to an equivalent argument in cycles, Y, and the 
OB3A 344 ; octant, I, contanining the argument. Y is then evaluated in a 
pB3A £9 3 polynomial chosen as a function of I to compute HSIN(X). 
OB3A re MTHSHSIN_RS:: 
50 73FD OB3A 70 TSTH RO : Check the sign of RO 
10 18 pB30 704 BGEQ POS_SIN : 
OOOOOB4A'EF 16 OB3F 705 JSB IN 3 RO/R1 = HSINCIX!) 
50 50 72FD OB45 706 MNEGH RO, RO ; RO/R1 = HSIN(X) 
05 0849 707 RSB 
Bit 8 
50 8000 8F AA OB4A 710 BICW #*x8000, RO 3; RO/R1 = ix! 
OB4F 711 POS_SIN: 
50 F4AC CF 71FD OB4F oe CMPH H_PI_OV_4, RO ; Compare pi/4 with iX! 
14 0855 71 BGTR SMALC_SIN ; No argument reduction is necessary 
50 F4B4 CF 71FD O0B57 714 CMPH H_9 PT_OV_4, RO ; Compare 9*pi/4 with IX: . 
3¢)60d«19'«OBSD) Ss 715 BLSS LARGE_SIN 3; Use special logic for iXi > 9pi/4 
at 
tH oe 3 pi/4 =< IX < 9epis4 
5E 20 C2 OBSF 720 SUBL #32, SP ; Allocate 32 Longwords on the stack 
$868 ts) ; for the reduced 
OOOOODA3'EF 16 O0B6 7 § JSB REDUCE _MEDIUM 3; Medium argument reduction routine 
0B68 7 : Y =_reduced argument on the stack, 
0B68 724 ; R2 = octant, R5 points to YHI and 
at 725 : _ R4& points to YLO 
07 O01 52 8F 086 ? § M_SIN: CASEB R2, #1, #7 ; Branch to one of four polynomial 
ee : evaluations depending on the 
6DE' ee 8 1$: «WORD P_COS_R-1$ 
6DE' 7 -WORD P_COS_R-1$ 
7A0° 0B 730 «WORD N_SIN_R-1$ 
7A0° 087 731 «WORD N_SIN_R-1$ 
73C* 0874 7 § «WORD N_COS_R-1$ 
073C° 1 7 - WO N_COS_R-1$ 
7AA’ 087 734 -WORD P_SIN_R-1$ F 
7AA® OB7A 735 ° P_SIN_R-1$ : octant bits. 
bare 739 
Oat 4 8 3 Logic for small arguments. iXi < pi/4. 
B7C 740 ° 
bare 741 SMALL_SIN: 
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E 
QO000E3D "EF 
54 

14 

07 00 52 
52.1 

FF 8F 

03 00 52 


; Floati 
ATHSHSIN_ 
81 7¢ 
18 OB 
7DFD 0B 
S4FD 
SFD 8 
64FD 0B 
OF B96 
S Be 
B98 
B9B 
C3 9B 
D BIE 
16 OBAI 
OBA? 
D5 OBA? 
13 A9 
AB 
8F OBAB 
085C° OBAr 
078D' 1 
078D' 0883 
tt 5 
0853° 7 
Se Se 
0853" OBBD 
ee 
ee 
8A OBBF 
9C pate 
8F OBC 
0BCB 
0987" OBCB 
0993" OBCD 
098C' OBCF 
0993° OBD1 
0BD3 
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DEGENERATE _CASE_ 


1$: 


oR 91:39:4 cr 
6-SEP=1984 11:25:20 CMTHRTL.SRCIMTHHSINCO.MAR; 1 


#*X3FC7, RO 
RO, -(SP? 
0: #SINLEN 


Reb, : 
(SP)+, RRo 


#32, SP 
REOUCE_LARG 


RG 
DEGENERATE_ 
R2, #0, #7 


ad ead cae a en ced ed od 
PRAAMAAMMM 


P_ONE-1$ 
UNFL -1$ 
N_ONE-1$ 
URFL -1$ 


R=-1, SINTBR 


Sete Ge Ge Ge Be Se 


CASE_SIN 


1 YAX/VMS Macro Vv04-00 


Compare with 2°-57 


We Wy mites evaluation is needed 


ao 
ow 
~v 
r] 
2. 
>< 
cm 


i = 
RO/R3 = 
= HSIN(X) 


attocete 8 tovby RS on the stack 


Be" nted 

; (R5) = YO.” (RE d= 
; R2 = octant bits 
; Check for degenerate case 


= YHI 


Compute index as (R2 - 1)/2 


| 
! 
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THSHSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:39:11 VAX/VMS Macro V04-00 Page 20 
-007 ATHSHSIN. RS or Ei obe 94:38:35 UMTMRTL. SREIMTHHSINCO.MAR: 1 ° (9) 
D 784 
BD 785 
+4 q § -SBTTL MTHSHCOS_RS 
BD 788 ; This routine computes the HCOS of the G-format value of RO/R1. The 
BD 789 ; penputes ten is performed one of three ways depending on the size of the 
44 £30 ; input argument, X. The processing is the same as described for MTHSHSIN_R4. 
D 036 
D 793 MTHSHCOS_RS:: 
50 73FD OBD 794 TSTH RO ; Check for reserved operand 
2p £008 SF AA D60«=—795 BICW #*x8000, RO ; RO/R1 = ix! 
0 F4 F 71FD DB £38 CMPH H_PI_OV_4, RO : Compare pi/4 with ix! 
5 14 —1 79 BGTR SMALC_COS 3; No argument reduction is necessary 
50 F428 CF 71FD OBES 798 CMPH H_9_PT_OV_4, RO : Compare 9*pi/4 with iX: 
5A 19 OBE9 799 BLSS LARGE _COS ; Use special logic for ix! > 9*pi/4 
OBEB go 
OBEB 01 ; 
OBEB 506 3 pi/& =< IX! < 9epi/4 
OBEB $ ; 
SE 20 C2 OBEB 04 SUBL #32, SP ; Allocate 8 Longwords on the stack 
OBEE 05 ; _ for the reduced argument 
OOOOODAS'EF 16 OBEE Be JSB REDUCE _MEDIUM 3; Medium argument reduction routine 
OBF4 0 : Y = reduced argument on the stack, 
OBF4 808 3 R2 = octant, R5 points to YHI 
OBF4 809 3 and R4 points to Y 
07 #+O* «52 8F OBFS 810 M_COS: CASEB R2, #1, #7 ; Branch to one of four polynomial 
OBF8 811 3 evaluations depending on the 
0714 OBF8 Big 1$: .WORD N_SIN_R-1$ 
0714" OBFA 81 [WORD NZSINR-1$ 
0680" OBFC 814 «WORD N_COS_R-1$ 
0680" OBFE 815 “WORD N7COS"R-1$ 
O71E' 8683 16 «WORD P_SIN_R-1$ 
O71E* OCO 1 «WORD P_SIN_R-1$ 
0652" 0C04 818 “WORD P7COS"R-1$ 
0652° 0C06 819 -WORD P_COS_R-1$ ; octant bits. 
0C08 0 
OC08 852 | Logic f LL arguments. :X! < pi/6 
3 Logic for small arguments. ‘Xi < pi/4. 
5608 g : 
coe 4 
CO 5 SMALL_COS: 
50 4000 8F 81 C08 § CMPW #*x4000, RO ; Compare 1/2 with {X: 
1 14 Bcoe BGTR 1$ : Sufficent overghang is available 
7E 50 7DFD OCOF 8 mOVO RO, -(SP) 3; (SP) = xX 
08 AE FFFFOIFF 8F CA 0C13 9 BICL @*XFFFFOIFF, 8(SP) ; 
OC AE 04 OQC1B 0 CLRL 12(SP) 3 (SP) = XHI 
7E 50 E 63FD 3 1 SUBH3 (SP), RO, -(SP) : (SP) = XLO 
4 3 dD0 OC : MOVL. SP, R4 : R4 = pointer to XLO 
oe 1 C26 ADDL3 #16, R4, RS ; RS = pointer to XHI 
Q6ga 31 OCeA 4 BRW NEEGS oust é ; Use special logic to obtain overhang 
50 3FC7 8F 81 C 3 5 18: CMPW #*x3FC7, R ; Compare with 2*- 
C 8 OC § BGEQ 2s 3 No polynon’ al evaluation is needed 
50 0 Gord C34 MULH2 aang 3; RO/R1 = X*x 
F561 CF OC 0 75FD 0C38 8 POLYH RO, #COSLENR1-1, COSTBR1; RO/R1 = HCOS(X) 
nh. * 


ee) 
2 


ee te, 
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MTHSHSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:39:11 VAX/VMS Macro v04-00 Page 21 MT! 
2-007 ATHSHCOS. RS o-$Ep=1 984 91:38:30 MTHRTL.SRCJMTHHSINCO.MAR;1 fay 2- 
50 08 70FD cag 41 28: MOVH #1.0, RO : RO/R1 = 1.0 = HCOS(X) 
05 0C4 4¢ RSB 
C4 4 
C4 44 
C4 45 LARGE_COS: 
5 ¢0 3 C4 $ SUBL #32, SP ; Allocate 8 Longwords on the stack 
Be C4 4 MOVL 3 ointed to by 
QOOODESD'EF 16 OC4B 848 JSB REOUCE_LARGE : (R5) = YLO, (R&) = YHI 
c51 49 ; R2 = octant bits 
54 rf bee 50 LCOS: TSTL RG ; Check for degenerate case 
“4 6} 51 BEQL § DEGENERATE_CASE_COS 
07 00 52. BF OC 2¢ CASEB R2, #0, #7 
OG * 0¢59 «854 18: .WORD P_COS_¢€-1$ 
B2' C38 55 .WORD P7SIN“C-1$ 
7A9* OC5D 26 -WORD N-SIN7C-1$ 
743" OC5F 5 WORD N-COS-C-1$ 
0743" 061 38 WORD N7COS"C-1$ 
07A9° C63 85 .WORD N7SIN7C-1$ 
07B2' 0C65 860 -WORD PUSIN7C=1$ 
06E3" 0C67 861 .WORD P~COS"C-1$ 
0C69 6¢ 
C69 86 
C69 864 DEGENERATE_CASE_COS: 
2 BA 0C69 866 BICB #1, R2 :; Compute index as (R2 - 1)/2 
52 52 FF BF 9C ogc 867 ROTL #1, R2, R2 
7 a 8F 0C71 868 SEB R2, #0, #3 
0C75 869 
O8E9" 0C75 870 1$: .WORD UNFL -1$ 
O8E2' OC77 = 871 .WORD N_ONE-1$ 
OBES" 0C79 872 -WORD UNFL -1$ 
O8DD' 0C7B 87 -WORD P_ONE-1$ 
0C7D 874 


palin saeabe 
| 


H 8 
ing Point Sine, Cosine and Sincos 16-SEP-1984 01:39:11 VAX/VMS Macro V04-00 Page MT! 
NCOSD_RS gr sera t9 4 94:38:3 MTHRTL.SRCIMTHHSINCO.MAR; 1 ° (98) é-' 
4 16 «SBTTL MTHSHSINCOSD_RS 
gp 4 ; This routine computes the HSIND and HCOSD of the H-format value of RO/R3. 
D 79 ; The computation is performed one of two ways depending on the size of the 
gp 9 ; input argument, X 
70 § 3 1) If iki < 465, then X is used directly in polynomial approximation 
70 : of HSIND and HCOSD. 
70 4; 2) If 45 =< ixt, then the subroutine REDUCE_DEGREES is called to reduce 
70 5 3 the argument to an equivalent argument in degrees, Y, and the 
70 § 3 octant, I, containing the argument. Y is then evaluated in two 
70 5 polynonials chosen as a function of I, to compute HSIND(X) and 
p : 3 COSD(X). 
7D 890 MTHSHSINCOSD_R7:: 
50 73FD 7D 91 TSTH RO 
10 18 3 BGEQ SINCOSD 
50 8000 8F AA 9 BICW #*x8000, RO 3 RO/R3 = ix! 
00000C92"EF 16 94 JSB SINCOSD 3 RO/R3 = HSIND(IX!) 
D 95 : R4/R7 = HCOSDC(IX!) 
50 50 72FD 8D 9 MNEGH RO, RO 3 RO/R4 = -HSIND(IX!) 
05 4 , 
9 99 SINCOSD: 
50 F3D9 cf 71FD 9 00 CMPH H_45, RO 3; Compare 45 to iX! 
0 14 9 901 BGTR SMALL. SINCOSD ; special processing for small arg 
5 10 C2 9A 90 SUBL #16, SP 3; Allocate 4 Longwords on stack 
5 SE D0 OC9D 90 MOVL SP, R 3; RS points to octaword on stack 
00001197°EF 16 0 a 2 : JSB REDUCE _DEGREES : a = he argument 
; = octan 
7E 5 7DFD A6 906 MOVO 5), ~(SP) 3; Save reduced arg 
0 BB OCAA 907 PUSHR #*M<R4, R5> : Save_octant bits and pointer 
00000078" EF 16 OCAC 908 SB EVAL_COSD : RO/R3 = HCOSD(Y) 
0 BA B2 909 POPR “m<R4, RS ; R4/R5 = octant bits/pointer 
65 6€ 7DFD O0CB4 910 MOVO (SP), (R5) : (R5) = octant/reduced argument 
6 50 7DFD OCB8 911 MOVO , (SP) : Save HCOSD(Y) 
00000014 "EF 16 BC aig $B EVAL_SIND 3 RO/R3 = HSIND(Y) 
54 7DFD C2 91 MOVO (SP)F, R4 3 R4/R7 = HCOSD(Y) 
4 co C6 914 ADDL 6, § : Pop 4 Longwords 
05 OCC 915 RSB 
Ocea 919 
CA 318 SMALL_SINCOSD: 
CA 91 UBL pe SP ; Allocate 8 longwords on stack 
CD 86920 mOovO RO, (SP) 3 Save_argument 
D1 0921 SB SMALL_COSD : RO/R3 = HCOSD(:X!) 
o7 «9 § MOVO RO, 186(SP) 3; Save HCOSD(:X!) 
o¢ 9 MOVO = (SP)+, RO : RO/R3 = argument 
E 924 $B SMALL_SIND ; RO/R3 = HSIND(X) 
E6 925 MOVO )+, RG 3; R4/R7 = HCOSD(IX!) 
EA 926 RSB 


ane i | 


Ne 04-00 Page 23 mT} 
; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:39:11 YAX/VMS Macro V wa: Pi] , 3 =( 
— ATHSHSIND RS gee b= o8e 94:38:30 MTHRTL.SRCIMTHHSINCO.MAR; 1 6 
cee 3 -SBTTL MTHSHSIND_RS 
ces 
; he HSIND of the H-format value of RO/R3. The 
EB ; ; cemalaitan te eerteresd ont of two ways depending on the size of the input | 
CEB § 3 argument, X 
EB 4 ; 1) If iki < 465, then X is used directly in polynomial approximation 
CEB «(935 ; of HSIND. é 
; ixt, then the subroutine REDUCE_DEGREES is called to reduce 
‘EB 3 $ : ” oy Ba to on equivalent argument in degrees, Y, and the 
cEB 3 ; octant, I, containing the argument. Y is then evaluated in two 
CEB 4 5 : polynomials chosen as a function of I, to compute HSIND(X). 
CEB 40 
50 73FD CEB We a : RO/R3 = X 
ee ae a rn et ad's dean 
36 § 72FD OCF6 945 MNEGH RO "RO ; RO/R3 = -HSIND(:X!) 
oft! EB: Reet 
CEB 969 on ercw —#*x8000, RO ; RO/R3 = ix! 
en 950 POS_SIND: 
500 951 7 CMPH H_45, RO ; Compare 45 to {X! 
vores © 95 BGTR SMALL_SIND 3 special processing for small arg 
0 OS bos ¢ SUBL #16, SP 5 Allocate 4 longwords on stack 
2 My 6 Op0B 934 MOVL SP, R 3 RS points to ar teuerd on stack 
00001197°EF 16 ODdE $36 JSB REDUCE _DEGREES ; cos? = reduce argume 
D14 95 
D14 958 EVAL_SIND: 
7 oats" D8 360 1$ COORD PSN b-18 
738° ODIA 961 “WORD P7COS"D-1$ 
$F * OD1C 996 -WORD P_COS_D-1$ 
0813° 0D1 368 -WORD P -SIN"D-1$ 
B0F* OD 964 .WORD N7SIN7D-1$ 
7A1° OD 965 -WORD N_COS_D-1$ 
7A1* 0024 96 -WORD N_COS_D-1$ 
80F* : 6 $68 - WOR N_SIN_D-1$ 
if 36 
p opy SMALL SIND: ; 180/pit2*-57 with ix: 
ol hich hig’ Ee 353 Bere events 1 ; cope tne evaluation is 
6 7rd OD 5 97 MOvVO RO, -(SP) : mecqueary 
g ¢ et D 4 974 MOVL SP, RS ; RS points to argument 
orf 73 H i $76 1$: FETA BSN. : Check for zero 
UF THF Sees, eAPh ie : Check ssible underflow on 
wae ee Oe oe Gree | Bg MLEST_PEG. 0; Chtonyersion fo rediens 
ost 31 bt? 3 9 BRW UNFL : Under tow will occur on conversion 
= . s 
so Piso oF Bip Bee i a Fe How Eh no ABMha pis ag0 - arabreia 
SUBW #*X6, (SP D4 . tA _ ; 2 rT 
ie 60rd be 3 4 ADDH2 (SP)+, RO : RO/RS = HSIND(iX!) = (pi/180) ix 


MTHSHSINCOS 
2-007 


Arhsnsi 


ing Point Sine, 
NDORS 


05 OD5B 985 3$: 


Cosine ond’ a ya 4 


RSB 


=1986 91:25:20 


AX/VMS 
MTHRTL. 


v04-00 


Pa 
SRTHHSINCO.MAR:1 72° 


Kk 8 
MTHSHSINCOS Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:39: AX/VMS Macro Vv04-00 Page 
2-007 ATHSHCOSD ORS pane 7 94:38:34 ERTHRTL. SREIATHNSINCO.MAR: 1 . dj ) 
4 f 7 -SBTTL MTHSHCOSD_RS 
p5¢ § 3; This routine computes the HCOSD of the H-format value of RO. The computation 
DSC 0 ; is performed one of two were depending on the size of the input argument, X: 
: 91 ; Details are given in the discussion on MTHSHCOSD_R4. 
psc 38 MTHSHCOSD_RS:: 
50 73FD ODSC 94 TSTH RO 3; Check for reserved operand 
50 £000 BF AA ODSF 995 BICW #*x8000, RO : RO/R3 = {Xi 
50 =F30 a vite +44 4 ast. Hrted hse ; Compare 45 to iX! 
C3 D6C 38 SUBL #16, SP ; Allocate 4 Longwords on stack 
SE 0D D6F MOVL SP, R : R points to octaword on stack 
00001197°EF 16 OD 1000 JSB REOUCE _DEGREES 3 (RS) = reduced argument 
4 ! 1 3; R4 = octant 
p78 1 § EVAL_COSD: 
07 OO 54 8F OD78 1004 CASEB R4, #0, #7 
06D7" OD7C 1005 1S: <WORD P_COS_6-1$ 
7AF* OD7E 1 $ «WORD P_SIN_D-1$ 
AB’ 0080 100 -WORD WN_SIN_D-1$ 
73D° O08 1008 -WORD N_COS-D-1$ 
73D" 0D84 1009 .WOR N_COS_D-1$ 
7AB' +t} 1010 .WO N_SIN_D-1$ 
O7AF* OD 1011 -WORD P_SIN_D-1$ 
06D7" OD8A 1 \§ «WORD P_COS_D-1$ 
BE ie 
D8C 1015 SMALL_COSD: 
50 F2FF CF 71FD $35 1oig CMPH H_SMALLD, RO 3; Compare 180/pi*2*-57 with ix: 
A 14 QD9 + 9 BGTR 1 3; Check if polyinomial evaluation is 
43 0 7DFD O0D94 1018 MOvO RO, -(SP) 3 necessary. 
5 E Bo D98 1019 MOVL SP, RS 3 RS points to argument 
068 1 1020 BRW P_COS_D 3; POLY needed 
50 08 70FD OD9E 1021 1$ MOVH #T, RO ; RO = 1, = HCOSD(IX!) 
05 ODA 10 ¢ RSB 
ODA 102 
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HCOSD_RS 66-SEP-1984 11 MTHRTL.SRCJMTHHSINCO.MAR; 1 
DA3 1025 
DAS 1099 
DAS 1 8 -SBTTL REDUCE_MEDIUM 
A 
DAS 1 ? ; This routine assumes that the absolute value of the argument, X, is in RO/R3 
DAS 1 § 3 and that pi/4 =< iX! < 9epi/4. It returns a pair of H-format values for the 
DAS 1 3 reduced argument in eight lLongwords on fhe stack previously allocate by the 
DAS 1034 ; the calling rout ing. hese longwords will be denoted by TO (4 gtf the stack 
DAS 1035 ; pointer? through 17 (32 off the stack pointer). YHI is in 14/T7, and YLO is 
A } § : in TO/T3. The octant bits in are returned in R2. 
DAS 1 3 : The reduced orgunent is obtained by locating the octant that X is in through 
DAS 1 3 a binary search and then subtracting off a suitable multiple of pi/2 
aA 18 | 
DAS 1 $s REDUCE MEDIUM: 
55 if 14 C1 QDAS 104 ADOL3 #20, SP, R5 3 RS pgints to 14 
65 50 7DFD ODA? 1044 MOVO RO, (RSS : 14/17 = X 
50 F280 CF “ep DAB 1045 CMPH = H_5_PI_OV_4, RO : 
pBi 1 r BLEQ 5$ 3 IX! >= Sepis4 
50 F208 SE BBE Tee BReLYPEOH AO sera 
52 01 00 et 1049 MOVL #1, R2 First quedrent 
17 «11 Ope 30 BRB SUBTRACT 
52 03 v0 ODC 1936 3$: MOVL #3, R2 3; Second quadrant 
12 #11 4: ' 27 BRB SUBTRACT 
so F276 ge rip BOIS 1055 58: CW WAT PLO MD lg 
52 05 0 Opcp 1987 MOVL #5, R2 : third quadrant 
0 11 DD 1058 BRB SUBTRACT 
52 07 00 oDp 1060 7$: MOVL #7, R2 :; Fourth quadrant 
00 11 ODDS 1061 BRB SUBTRACT 
By 1068 
DD7 1084 SUBTRACT: 
54 04 1 OQDD7 1065 DL3 #4 4 R4 3; R4 = points to T0 
53 FD Az4g 3E DDB 1 66 MOVAW =-3(R Cred R3 ; R3 = index into PI_OV_2 table 
53 F2kB CF4 E ODEO 106 MOVAQ PI_OV_2CR3i, R3 : R3 =_pointer into PI_OV_2 table 
6 03 6erp 434 ! $8 susie (R3)+> (RS) : 14/17 = Yar appreniads ten to YHI 
D6 ODEC 1 19 INCL. R : Adjust octant bits 
3.6COS $000 or AB OQDEE 1071 1$: BIicw3 9 #*x 000. (R5), R1 3; R1 = exponents bits of ‘YHI'! 
51 F95 BF 6©B1 ODFS 1 i CMPW #*XxX3F95, R1 : Check for at least 6 significant bits 
18 14 +h ! Y BGTR NOT_ENOUGH_BITS 3 in YHI* 
64 $3 7D0FD ODFB 1075 MOoVO (R5), (R4) 3; T0/T3 = YHI' 
08 AS FFFFOIFF BF CA ODFF 1 id: BICL #*XFEFFOIFF, 8(R5) ; 
OC A D4 OE07 107 CLRL 12(R5) : 14/17 = high 56 bits of YHI* = YHI 
64 6 ogre EQOA 1 A SUBH Hg (R4) ; TO/T3 = low bits of YHI' 
635 6 Fe E 107 SUBH (R35), (R4) 3; R4/R5 = YLO 
3 Be IRS 


MTHSHSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:39:11 VAX/VMS Macro V04-00 Page 27 MT 
2-007 REDUCE.MED 6- Pa} 9 33:3) MTHRTL.SRCJMTHHSINCO.MAR; 1 . (79) é- 
? eB. 5 NOT_ENOUGH_BITS: 3; Add an apert tenet 49 significant bits 
€1 1 3 to YHI ° att YH 
7E $3 7DFD OE1 1 nove (R3), =(SP) 3: (SP) = 113° additional bits (H-format) 
08 AE 7C O17 1 Q CLRQ. = (SPS : (SP) = 49 additional bits (H-format) 
65 oF 6 FD OEFIA 1 : SUBH2 (SP), (R5) : 14/17 = 
64 ES FD OEE 1 SUBHS }§ (SP)#, (R3)+, (R4) : T0/T3 = oe 64 bits of YLO 
6 63 6 FD OE23 1 8 ADDH2 (R3), (R4) : T0/T3 = -YLO 
64 8000 at Eo? 1 XOR #*x8600, (R4) : T0/T3 = YLO 
05 Ofc 1 : RSB 
E2D 109 
E2D 1092 
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0000003D 
00000020 
00004000 
0000 7 

00003F 


00000000 00000000 00000000 00004011 


54 50 OOOOFFFF 8F 
54 (01 
OC AS 8654 = «610 


CB 
8B 
9C 

06 AS 51 10 9¢ 
03 1 

0c A506 

O46 AS 52109 
03 1 

08 A506 

65 53 10 9 
03 1 

04 a5 06 

50 3FF2 BF A2 

9C 

st FRPFFOOS BF CA 
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om 
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oo 
mo 
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This routine is used to reduce vere arguments (iX! >= 9*pi/4) modulo pi/4. 
It returns the reduced argument, Y, in R4/R5 in units of cycles, and returns 
: the octant bits, I, in R2. 


+ The method of reduction is as follows: 


x#(4/pi) = anetsta(et where nis an integer and 1/2 =< f < 1 

; = 2*(m-113) #(2*1158f) #(4/pi) 

: = (2°113%f) #( *insh13)84/p 

; = K*C, where K = “6 13*f 7 } - ae er yes = 2e(n-113)*4/pi 
: Let L = KeC estate 8, where 0 =< L < et he integer(L) and 
; h = fract(L), then if I is even Y = h, icles = 1-h 


; CONSTANTS: 


L-INT_ WEIGHT = *X3D 3; weights exponent by $3 
u “TERA_WEIGHT = “XxX 3 weights exponent by 32 
WIMAX_GEIGHT = *x4000 ; maximum unbiased exponent 
WIYLO-WEIGHT = *X4 ; petents gupeness | by 
W-ADJOST = “X3FF2 ; Used to locate b binary spoint in 
s MTHSAL_4_OV_PI_V tab 
H_2_T0_32: 
-OCTA “%X4011 s 2°32 
REDUCE LARGE: 
3; The rare step 4 to 9 gbtein the location of the o Btnery | point in the represen- 
3 tation of C = 2*(n-115)*(4/pi) in two parts - er of longwords from 
3; the start and the. mt ate of bits from the most "slanttiosal bit of the next 
3; longword. Also K = 2°1135*f must be obtained. 


BICL3 #*XFFFF, RO, R4 ; R4 = Ist 16 fraction bits of Xx 
BISB #1, R4 3; Restore hidden bit 
ROTL #16, R4, 12(R5) > 12(R5) = High 17 bits of kK 
ROTL #16, R1, 8(R5) ; 8(R5) = Next 32 bits of K 
BGEQ 1$ :; Check for high bit set and 
INCL 12(R5) s adjust if necessary 

1$ ROTL #16, R2, 4(R5) : 4(R5) = Next 32 bits of K 
BGEQ Hi ; Check for high bit set and 
INCL (R5) § adjust if necessary 

2s ROTL #16, R3, (RS) 3; (RS) = Next 32 bits of K 
BGEQ 3$ ; Check for high bit set and 
INCL 4(R5) : adjust if necessary 


3$ SUBW #W_ADJUST, RO ; Unbias om and adjust for leading 
: poreet, RO = location of binary 
ROTL #-3, RO RB : Divide. Ro by 32 and mull by 4 to 
BICL #°xXFFFFFOO3, R3 ; R3 = # of Lonquords (in bytes) 4 
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address of RTL vector entry 
eddress 9 MTH$AL_4_OV_PI table 
oints to 1st quadword of interest 
:0) = # of bits within. lLongword 


MOVAL G*MTHSAL_4_OV_ 
GATHSAL- 470 


R 
W*KE0, RO 
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The next step is to generate an aporen teat ten to 


C it C’* to be used 
e Comput ing wore will consist a the firs i 
r ts 


e integer bits of 
C and the first 6 heastten bits of C. These bit be obtained tren 6 


t 
constant stored in the interger array MTHSAL 4_OV_P 


NOTE: The ASHQ, ADDL, and MULL instructions in the follow sections may 
result in an integer overflow trap. The overflow incurred is intentional, 
; so that the IV bit must be turned off. The IV bit is not restored until 
; after all of the necessary fraction bits hav been generated. 


MOVPSL =(SP) ; Put current PSL on stack 
ICL #*C<PSLSM_IV>, (SP) ; (SP) = current IV bit 
BICPSW #PSLSM_IV ; Clear integer overflow bit 


mE 


7E 
6E FFFFFFDF 8F 
20 


Doo 
oro 


The necessary calculation to produce the reduced argument can requires up to 
26 longwords of oh eper ety, work space. This work space will be allocated 

on the stack. e 26 Longuords gh teested. only 16 are needed for the 
duration of the reduction scheme. for the Gyrpesse of comments ae 16 
locations will be refered to as TO through TE and will be accessed by 
pointers in R R4, and SP. The stack will be allocated (after initial- 
ization) as folteut: 


SHHHOGSOOGSO OS OSSD DOS OVOSOOSOSOSOOOOONONa BD 


PDMOOOOO OOOO OO 060 69 09 09 09 09 09 09 08 NINN NNN SSP PAA AA AAAA MIA 
NOAM WN @ SO ODONAU EWN 9 OD NAUE WIN (0 ODNA UE WH OUOONOUS 


Pp 
E 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
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: T0: 
: Temporar _storage for use } 11: | 
3 Sucig ¢ eneration of : | succcccascencoesonss ' 
3; addit ~ J bi s of the pro- | T2: | H 
; duct KeC ' | eceenscenesseneenses ‘ 
; | 13: | | 
4 : T4: 
38 ; Additional bits of c’ =| ‘TS: | 
38 ; t 16: | : 
98 ; t 17: | 
Be OS pee | 
98 3 : Contains all of the ge ; ; 
9B 0; significant bits . ‘ ‘ ° 
98 1; generated through . ‘ . . 
98 § 3 current iteration H 
98 1208 1D: | ! 
98 1209 Te: | 
98 i : Weighting factor 1 TF: } 
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5E 00000068 8F 
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mMOVa RO, (SP) ; Save RO/R1 
; Set up pointers for multiplying 

ADOL #36, sf. R1 ; R1 points to Ist longword of C*' 

ADDL #32, R2 : R2 points to T7 

ADDL #24, RD RG ; R4 points to TD 

CLRL ne} This code generates_the product 

EMUL (R5) we! | Yi (R2), (R2) of the, iow order 32 bits o 
CBL and stores them in 


oure> ff 
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Get C'* = €(0):C(1):...:C(7) on stack. C(0) though (C€(7) 7pre a tenes 

sateners generated from the bag 8 representation of C. The high three bits 
of C(O) are he. the first three bits to the left of the bin 53 Bist of C 

The remaining bits C(0) and CC) ehouss C(7) are the first its to as 
right of the binary point of C. Note that the C(i)'s are adju 

ompensere for their signed (rather than unsigned) interpretation, in mts EMUL 
nstruction. Note also that the representation of C has no more 

consequtive ones, so that no carry is possible from the aistesele 


SUBL #104, SP ; Allocate 26 longwords on the stack 
ADDL3 o36° SP, Re ; R2 points to longword before C(7) 
ADDL3 #28, R2, R ; sntters tee one S counter. R3 points 
; _ to longword before 
ASHQ RO, (R1), (R3) : add + Sener proper F quedvord so that 
SUBL #4, R3 ; R3 a ll to te word before ((6) 
SUBL #4, R1 ; R2 points to next quadword in 
; THSAL_4_OV_PI_V table 
ASHQ RO, (R1), (R3) ; Shift quadword so “that C(n) is in 
; the sppropr late Longowrd 
BGEQ 5$ ; Check for high bit of C(n) set 
INCL 8(R3) ; Bit set. Adfust t etn -1) 
ACBL R2, #-4, R3, 4$ 5 Loop umest 0) though C(7) are on 
3 e stac 


Generate the low 256 bits of the product K*C'' = L. This eregucs is 
equivalent to multiplying K times C'* modulo 8. The result of t 
cgduct is in T7/TE with bits 31:29 of TE the octant bits, and the remaining 
bits the faction bits of the product. The last 113 fraction bits (bits 
os non-valid fraction bits that will be used later if more fraction bits 
need to be generated. 


A 4 6$ c'* ang 's 
MULLS (R5) tris’ 4(R2) location 17/TE 


DDL  4(R25, (R25 


CLRL R3 ; R3 is the loop counter 


| 
| 
eet LC CC CC LL LL 


aetna mamma mamma amma aameaaamaaataaammaamaaamaaaamaaaaaaaaamaaaaaaaataaasaaaaaaaamaaaaaamaasaeaamaa eames acamamaamaaacamaaasaaamaiamaiaatag a | | real 
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EEE 1265 
2 98 AE43 bf EEE 1 06 7$: MOVA ocsr)tnss RO ; RO and R4 point to Ist and Last 
4 5E iC C1 OFFS 126 ADDL3 #28, SP, R4 :  lLongword of temporary storage 
EF? 1 rt: : for computing the product of the 
EF? 1 § 3; _.next 32 bits of | and C'' 
22 % CO OEF? 1270 ADDL #4, RS : RS points to next 32 bits of K 
51 5& 2 C1 EFA 4 ADDL3 #36, SP, R1 ; Ri points to 1st Longwork of C'’ 
$0 pf EFE 1 ie CLRL Hf 3; This loop oggerotes the product of 
4 3 A 4. 1274 8$: EMUL (R5), (R1)+, (RO), (RO) ; the next 32 bits of K and C''. 
FFF5 50 04 4 Fi OF 1275 ACBL R4, #4, RO, ; The result is stored n+2 longwords 
AE 65 ¢5 FOB 1 06 MULLS (R5), (R1), 100(SP) : off SP, where n is stored in R3 
60 64 AE (CO F10 , ADDL 100<§P), (RO) ; 
50 08 AES DE OF 4 1 rg MOVAL oc seacRad RO ; RO points n+2 Longuors off SP 
52 4B AEG ms F19 1280 MOVAL 72(SP)CR34, R2 > R2 points to T7 + n Longwords 
54 04 C Ore ! + ADDL #4, R ; 
82 80 CO ors) 1 Hy ADOL (RO)+, (R2)+ ; This loop adds the product bits of 
82 80 08 OF24 1284 9S: ADwC (RO)+, (R2)+ 3 the previous loop to the bits 
pest 1285 3 glready generated. The results 
FFF7 50 00 54 F1 OF27 1286 ACBL R4, #0, RO, 9$ ; This curious construction is used to 
prs } fi 3; avoid effecting the carry bit 
FFBB 53 O01 O02 FI 8 ! +4 ACBL #2, #1, R3, 7$ 3; Increment loop counter and branch 
50 7D OF33 1291 mova (SP), RO ; Restore RO/R1 
55 OC C2 OF36 1636 SUBL #12, R5 ; RS points to K 
5E OF39 129 ADOL #40, SP : Pop 10 longwords off stack 
aia 
OF3C 1296 ; At this potas there may or may not be enough valid bits in R3/R4 to generate 
OF3C 1297 ; Y. If the first 12 fratction bits are all 1's or O's, there a possibility of 
OF3C 1298 ; loss of significance when coneyt ing Y.. Consequently, we must check for loss 
OF3C 1299 ; of significance before converting 14/T7 to Y and I. 
OF3C 1300 ; 
OF ee 4 
OF : : 8 3; Set up pointers 
53 SE 1C C1 et 1304 ADOL #28, SP, R3 ; R3 points to T7 
54 SE SC (C1 eee : 05 ADDL #60, SP, R4 ; R4 points to TF 
64 FC A4 0000080 8F C1 OF44 1 09 ADOL3 #*X80, -4(R4), (R4) : If the first 22 fraction bits are 1's 
64  3FFFFFOO a 03 F4D 1 08 BITL  #*X3FFFFFOO, (R4) : and the reduced arg = 1-f or the 
, hCG F541 BNEQ CONVERT ; first 22 bit are 0 and the reduced 
OF56 1310 : arg = f, then (and only then) bits 
F56 «(1311 : 29:8 are 0 and significance will 
F561 \§ : be lost. 
F56 «131 
ee 1314 ; 
36 1 1g 3; More bits neea to be generated to cover the loss of significance. 
F56 «1317 * 
OOO0108F'EF 16 OFS6 1 + JSB GEN_MORE_BITS ; Generate 85 additional bits and add 
F5C (131 ; them to existing bits. Results are 
= 3 $9 ; . stored _ in ; 
53 046 c2 OFS5C 1 SUBL #4, R3 3; Adjust R3 to reflect the addition of 


—————ee —y 
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FSF 1 ¢ ; another Llongword of K*C'' 
15 FC AG) «61D )=«6©EO «COOFSF 1 BBS #29, -4(R4), 108 ; Check if loss of significance is d 
+71 ! : ; to leading ones or zeros 
foe : § ; Lost significance due to leading zeros 
64 FC AS 08 00 FA F641 8 FFS #0, #8, =4(R4), (R4) ; If at least one bit is set. This | 
1 2 OF6A 1 BNEQ CONVERT 3; _.means lost a A ance was minor. 
FB AG OOOIFFFF BF 01 4) 1330 CMPL @*X1FFFF, -8(R4) 3; If one of the high 15 bits is set, 
17 15 OF74 1331 BLEQ co ; lost significance was minor. 
0084 «31 Hy ! § BRw LEADING_ZEROS 
ore : : ; Lost significance due to leading ones 
64 FC AG 08 00 £8 OF79 1336 108: FFC #0, #8, -4(R4), (R4) ; If at least one bit is clear. This 
Sm OFre Saar BNEQ CONVE ; _ meanslost significance was minor. 
FB AG FFFEOOOO 8F Di OF81 1 3 CMPL #*XFFFEQOOOO, -8(R4) :; If one of the three high bits is 
2 1& OF89 1339 BGEQU ONVERT ; clear, lost significance was minor. 
oc . a 4 4 1340 BRB LEADING_ONES 
FBD 3c 
OF 8D : “8 CONVERT: 
OF8D 1345 ; Isolate octant bits and convert fraction bits to a pair of D-format 
OreD 1347 quantities YHI and YLO 
64 FC AG =603.—s1D~Sté@EF OF 1348 EXTZV #29, #3, -4(R4), (R4) 3; TF = octant bits 
FC AS £0000000 8F CA OF93 1349 IcL_  #*x€0006000, -4(R4) ; Clear octant bits 
C3 OF9B 1350 SUBL3 #20, R4, R3 3 R3 points to low order bits of h 
QOOO1IOC'EF 16 44 132) JSB CVT_TO_H ; ites “Brake 
: = es 
10 AS BS OFAS 138 TSTW 16(R5) : 16(R5) could be zero. If so 
04 13 OFA8 1323 BEQL H_HI_0 3; .do not do SUBW #*X1D,16(R5) 
10 AS 1D Ae ae ' 22 ame SUBW #*X1B, 16(R5) 3 16(R5) = h_hi 
65 1D A2 OFAE 1357 ~ ~  SUBW #*X1D, (RS) : (R5) = h_lo 
62 . ore ! 38 by 189) ; Check for h_lo = 0 
0B 64 «EY OFB3 : 9 1$: BLBC (R4), 2% 3; Check for odd or even octant bits 
He : 66 3 Octant bits are odd. Reduced argument equals 1 - h. 
BB 1364 SUBH3) 16(R5), #1, 16(R5) : 16(R5) = YHI = 1 = h_hi 
FB 69 MNEGH (R5), (RS) ; (R5) = YLO = -h_lo 
He ' e 3 Get octant bits 
| 
FCS 1 $3 2$: MOVL (R4), R2 ; R2 = octant bits | 
FCG 1370 BRW GET_YHI_YLO | 
FC9 «1371 
Fea 1398 | 
FC9 1374 : At this point it has been determined that there is a major loss of 
FC9 1375 ; significance and the processing begins a looping phase. Each iteration of | 
FC9 1376 ; the loop will generate additional extra bits of K*C* until enough significant 
FC9 1377 ; bits to compute Y are available. Gur ing this time the nine longwords 
FC9 1378 ; allocated on the stack will be used as follows: 
| 
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TO/T4 Temporary storage used when generating extra bits. | 

T6/TE Contains all significant bits generated so far. | 

TF Contains a counter, W, indicating the appropriate exponent 

oe oe last lLongword of fraction bits used in converting | 

EADING_ONES: 


: at S Greceze=ne continues here it is known that the loss of significance is due 
a string of leading ones. 


MOVL WL_LINT_WEIGHT, (R4) ; TE = omp bias for soef Longword 
3 of the product KC 


; seen for . on significant bits 
; Enough bits. anes to floating. 
contains Ke 


LOOP_1: CMPL #*XFFFF8000, -8(R4) 
BGTRU CONVERT 1 
JSB GEN_MORE BITS 


Sete Ge Ge Ge & 
“ 
wm 
~ 
_" 


CMPL #-1> -B(R4) ; Check for all 1's 

BGTRU CONVERT. 1 ; Not all 1's. Enough precision bits 
; to compute Y 

MOVO =24(R4), a > Hrd 3 Compress representation 

MOVO -40(R4), Kec’ 


ACBW = #W_MAX_WEIGHT, #W_TERM “WEIGHTS (R4), LOOP_1 
: Increment “weighting factor. If 
3 ve ssgnting factor is greater than 
; then no more bits need to be 
3 Lowadh § 
The weighting factor is greeter than 16384. This means that the reduced 
argument is either not q {es tagusenenie from 1 or too small to be represented 
in F-format (i.e. underflow.) Zero is returned in R4 for the reduced 
argument to signal this occurance. Note that under these conditions the 
correct function value is one of the values 0, +/-1. The 
correc’ soetes is determined by the calling program based on the octant bits 
returned in R1. 


CLRL ; Reduced argument is zero 
aire 129, #3 - ~4(R4), R2 ; R2 = octant bits 
CONVERT_1: : 
SUBL3 #24, R4, R3 3; R3 points to low bits of h 
JSB CVT_TO_H 3 (R5)= 2*Wth Lo 
> 16(R5) = 2°Oeh_hi 
SUBH3 16(R5), H_2_T0_32, 16(R55 
; 16(R3) = 2°W*(1 = h_hi) 
» 16¢(R5) : 160R5) = 1 -h 
(R5) ; (R5) = -2*W*th_To 
3 check for h_lo = 
: (R5) = = h_To 
: R2 = octanf bits 


(2 


tT << 
m0 
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GEN_MORE BITS: 


H 
H 
H 
H 
H 
H 
H 
L 


0 

1 

: : This subroutine generates 145 extra fraction bits and adds them to the 

4 ; existing bits. : This routine is always entered via a JSB spetrys tten. 
5 ; Consequently, SP points to the first longword BEFORE T0, rather than T 

§ : itself. R3 points to the Last longword of previously generated bits. 

0 
1 
2 


ADDL3 #24, SP, R2 ; R2 points to storage for C(n) 

SUBL #4, RI ; Adjust potater to get next quadword 
: rom MTH$AL_4_OV_PI_V 

ASHQ RO, (R1), -4(R2) ; T5 = C(n) 


= 8&w & 
51 04 (C2 
FC A2 61 50 79 


atin 
E_LARGE 
1933 it 
! ° 1? 8 LEADING_ZEROS: 
102D 1440 ; If processing continues here it is known that the loss of significance is due 
102D 1441 ; to 4 string Hf eading zeros. Note that fe is known that the loop for 
102D 1408 ; leading zeros will terminate before an underflow condition occurs so that the 
! ° 1227 ; loop does not include a test for underflow. 
64 3d 00 19 D 1445 MOVL #L_INT_WEIGHT, (R4) 3 TE @ oxp bias for last longword 
: 1328 ; of the product K*C' 
FB AG QOOOIFFF BF D1 1 1703 LOOP_O: CMPL #*XOOOO1FFF, -8(R4) ; Check enough fraction bits 
~.) we .% 144 BLSS 3 Gough bits. Convert to floating 
QOOO108F*EF 16 103A 1450 JSB 1 MORE_BITS : T2/R7 contain KeC'' 
FB AG DS 1040 1451 18: TSTL -8(R4) ; Check for all 0's 
2 % + 1226 BNEQ covert .0 ; Not all 0's. Enough precision bits. 
EC AG =6EB AG 70FD «261045 =«1145 MOVO =24(R4), -20(R4) ; Compress representation 
DC AS DB AG 7DFD «4610468 «41454 MOVO =40(R4), -56(R4) 5 of KeC'' 
20 AO 1051 1455 ADDW #W_TERM_WEIGH : Increment weighting factor. 
DA 11 13ee 1528 BRB LobP_0 
54 18 (C3 1986 1289 CONVERT OeL3 #24, R4, RB R3_ poi low bi fh 
o RG, 3 oints to low bits o 
OOOO110C'EF 16 1928 1299 SB CVT_TO_H F (RSS oyeosh ns 
‘ = 2*We 
52 FC AG 03 1D EF 1060 1298 EXTZV #29, #3, -4(R4), R2 ; R2 = octant bits 
10 A564 C2 1066 146 SUBL  (R4S, 16(R5) : 16(R5) = h hi 
65 B5 106A 1464 TSTW (R5) ; Check for R_lo = 0 
03 13 106C 1465 BEQL GET_YHI_YLO : 
65 64 C2 106E 146 SUBL (R4T, (RS) 3 (RS) = h_lo HL! 
1071 146 H 
1071 1468 GET_YHI_YLO: HL 
54 55 10 C1 1071 1469 ADDL3 #16, R5, R4 ; R4 points to YHI H_. 
7E_ 64 7DFD 1075 1470 OVO (Ras, =(SP) ; (SP) = high bits of Y H_é 
ge He cate H8re Hi Sit EAE con spe se Sencar ee te 
° : a 
6 8E 60FD 94 1478 onstane ADDH (SP)+, (R5) ; (R4) = YLO He 
40 AE 3 1084 1475 BISPSW 64(SP) ; Restore IV bit and exit cal 
5E 00000044 8F CO 1087 1o08 DOL #68, SP ; Remove mask and temporary storage .! 
05 108€ 147 $B .! 
108F 147 .! 
108F 147 y 
108F 14 
108F 14 
108F 14 
108F 14 
108F 14 
108F 14 
108F 14 
108F 14 
108F 14 
1 ; 14 
1093 14 
1096 149 
1096 149 


H 9 
MTHSHSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:39:11 VAX/VMS Macro Vv04-00 Pa 5 
2-007 REDUCE. LARGE beeen 13ke $1:3ei20 EATURTS. SaeaarHingioco.maR:1 2% (30) 
33.18 } 38 1333 BGEQ 1$ ; Branch if high bit is clear 
! $0 1232 ; Logic to process unsigned values greater than 2°31 - 1 
ee ee ee en ee re 
OB AE OB AE, SE RE AE PA RA? 65h ERUL URS. 1H), cso, cso 
OC AE OC AE 6 85 é A A 84 13 1 ERUL ae , (R2), 12(SP), 12(SP) 
10 AE 10 AE 85 CA A 10C1 1 i EMUL (R55, RS) 16(SP), 16(SP) 
14 AE OCA 0 10C9 1504 ADDL } (R5), 20(SP) 
1E 1 4 : 5 BSB 
130 1307 3; Logic to process unsigned values less than 2°31 
04 AE O00 62 65 #7A 10D 129) 1$: EMUL (R5), (R2) "0. 4(SP) 
O08 AE OB AE 62 0405 7A 1006 1510 EMUL 4(R5S, (R25, B(SP), BSP) 
OC AE OC AE 6 08 AS 7A 10DE ae MUL 8(R5), (R2) 12(SP5 12(SP) 
10 AE 10 AE 62 OC AS 7A 1066 1 1 MUL =-:12(R5$, (R25, 16(SP5, 16<SP) 
aH 
+43 131? 3; Add new bits to old 
a2 S&S GS Ct 10F§ 1317 2$: ADDL3 #4, SP, R2 3 R2 points by tow order longword of 
3 of new s 
FC A ge DO 10F2 1519 MOVL (R2)+, =4(R3) 3; Move low order bits to end of List 
8 CO 10F6 1520 ADDL (R2)+, (R3) ; Add in new bits 
O04 A 82 pe 10F9 1521 ADWC (R2)+, 4(R3) : 
08 A HE D8 10FD 15 ¢ ADWC (R2)+, 8(R3) 3 
OC A 8 D8 1101 15 ADWC (R2)+, 12(R3) : 
04 1€ 1105 1524 3$: BCC 4$ 3; Check for carry from previous add 
82 D6 1107 1525 INCL (R2)+ :; Propagate carry 
FA 11 1109 15 g BRB 3$ 3 
05 1108 1527 4$: RSB 3 
110C 1528 
110C 1529 
ge as 
110¢ 15 é CVT_TO_H: 
Hare 
1106 1535 : This routine converts an array of five lLongword pointed to by R3 to a pair 
1198 1 § 3; of H-format values. The results are returned on the stack and are pointed 
110C 1 3 to by R5 
110C 1 8 
He 13 
3s HH SS 1149 1361 ADDL3S #12, R3, R2 : R2 point te second most significant 
50 SE 04 C1 EB ! ‘8 ADOL3 #4, SP, RO : RO points to T0 
6¢ 05 1114 1545 1$ TSTL (R2) ; Adjust for signed 
18 1118 1 6 BGEQ cb 3; =6 rather than unsigned 
04 4 b6 111 154 INCL 4(R2) : conversion 
FFEF 52 FFFFFFFC 8F 5 F 1438 ' $8 2s ACBL R3, #-4, R2, 1$ ; 


xz 
mm 
c-_ 
oo 
>s 
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ONN "PT PnP TPN 


0 CVTLH (R3)#, (RO) 
; SuBH AW TERM_WELGHT, (RO) 

3$: CVTLH (R3)4 
A Py." (ROD 
TERM_WEIGHT, (RO) 
d+, TSP) 

), (RO), (RS) 

BEQL «=s5$ 
SUBW © #W_YLO_WEIGHT, (RS) 
(R3)+, (RO) 
6 


SUBW HW Sa WEIGHT, (RO) 
oe +, TSP) 


Convert low order 
three longword to 
H-format and store 
on stack pointed to 
by R5 


60 


Pa 
Ld 
oww 
<cm 
_~ 
_ 
= 
Cota trian! 
ADEMM 
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65 0040 


Convert high order 
two lLongwords to 
H-format and store 
16 off R5 on the 
stack 


a 
A 
o 
a< 
m— 
= 
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2.10_32, (SP) 
tsby+. “(RO), 16(R5) 
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-SBTTL REDUCE_DEGREES 
; a routine assumes that the sbeotute value of the argument is in RO/R3. 


| 
; reduction process is performed in two grages. The first stage of 
; the reduction reduces the argument modulo 360 to a value less that 2°*112, 


a ed a a tt be a a a 2 a 2 tt td it mo 


DRL APAP DAA AA AAPA DDD DADA PADD Paras «= FTO 


ro V04-00 Page 3 
RC JMTHHSINCO.MAR; 1 (21 


the second stage reduces the argument modulo 45 to a value less than 45. 
; The reduced argument is returned in four Longuords previously allocated by 
; the calling program and pointed to by R5. The octant bits are returned in 


OTR rTe Tere rere 
@ ovr 
J 
a 


3; Constants used in this reduction: 


POWER_MOD_360_0: ; Powers of 2 modulo 360 for ti = 0 
0 001 -QORD 1, > 4, 
008 3938 4 010 «WORD 16, H ° 64, 12 
0130 00 100 -WORD 256, 152, 304, 248 
POWER_MOD_360_1: ; Powers gf 2 modulo 360 for ti <> 0 
8 0088 0110 0088 -QORD™ 136, 272, 184, 
0080 0040 it 4 0010 «WORD 16, : ° Gi. 12 
OOF8 0130 0098 0100 -WORD 256, 152, 304, 24 
REDUCE _DEGREES: 
50 4071 8F 861 CMPW #*x4071, RO ; Compare ix! with 2°112 
14 BGTR LAST_STEP 3; Branch to special logic for med arg 


It is assumed here that the argument is greater than 2°113. 
The argument is reduced as follows: 


BSDOOGOOCOOCOO OOOO OOOOOOWWDOOODWOODoo VII INN Wo 


CONAN WIN OS OD NAUES WIN OS OO ™! PUL WIN SO OD NAN EWN OOD NAN E WI  O ODNOAUE UM 


FF OOOO OOO OOD 000000 OO OO ODD SII INP AAAA AAPA AAA AAAS 
QeMM MMMM MMMMM MMMM AO VINNIE NNN NNN NNN 


meet ee et ts et a a sb — bs 2 2 a ss a ts a te mo 
me a a a at a os = a ss a 


$ be x = 2“tef, where t > 113 and 1/2 =< f < 1. And let J = 2°113*f = 
F *60*J1 + J2 and K = 2°(t-113). Since 2°60 = ontg modulo 360, we have 
3 J = 2°12*J1 + J2 modulo $0. Now Let 5, = t - 115 = 12*t1 + t2. Note 
3 that (2°12)%2 = (2*9)#(2"1 ) = (2*9)#(2*3) = 2°12 sotyie 360. Hence, if 
3 1 is not zero, K = 2°t' = 2°(12et14t2) = (2%12)#(2*t2) = 1 ane te ulo 
10; 60. Forti = 0 K = 2*t2 Sengoquenst . define Kk‘ congruent to 2*t2 if 
7 3 tl = 0 and congruent to 136*2*t2 otherwise, where 0 =< k' < $0. Then x‘ 
' : = K'#(2*12eJ1 + J2) is congruent to s modulo 360 and x* < 2113. 
54 50 00 1 MOVL RO, R4 3; R4 = exponent bits of X 
; 7FFFE 8F 8 =6AA 1 BICW @*X7FFF, RO ; Clear exp bits of X 
0 4071 8F A 1 BISW #*x4071, RO ; RO/R3 = J 
54 ee > | : SUBL RO, R4 : R4 = t’ 
63 0 7DFD A 1 OVO RO, (RS) : (RS) = J 
52 FFFFOFFF 8F CA B BICL @*XFFFFOFFF, R2 F 
3 D4 8 CLRL 3 ; RO/R3 = J1*2*°60 
5 62F BB SUBH2 RO, (R5) : (RS) = J 
0 2 A BF SUBW #*x30, RO 3; RO/R3 = 21505) 
50 65 60FD ce ADDH2 (R5), RO : RO/RS = 2°12eJ1 + J2 = J modulo 45 
65 54 C A? 3 DivwS #12, R4, (RS) 3; (RS) = tl 
—E 13 CA BEQL 1$ 
65 C A4 cc MULW #12, (RS) ; (RS) = 12*t1 


65 


65 


50 
65 50 


7E 3854 


54 65 


C16C 8F 
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65 
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eo 
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EP=1984 MTMRTL. SRCIMTHHSINCO.MAR: 1 Tat 
SUBW 3s (RS), RG . t2 
EVTWH POWER_ROD.360.1CR4], (Rs) 
BRB 2$ : (RS) = K 
1$: CVTWH POWER_MOD_360_0[R4], (RS) foal 
2$: MULH2 (R5), RO ? ROJR3 = Rae (mod 45) 0 =< RO < 2°53 
LAST_STEP: 


; Argument reduction scheme for arguments with absolute value less than 2°112 
The reduced orqueett Y is commuted as follows: 
Let 1 = {ne x/45) 
Ss even 


then Y = X = 45* 
else Y = (141)#4 


Be Gee Se Se Ge te ee 


3; Compare 2°36 with iX: 


; (RS) = iKi/45 


CMPW My ot RO 
Q NO_OVERF pow 


BGE 
MULH3 H_T_OV_4 (R5) 


Turn off IV to avoid an exception in EMODH 


are of RG 3; R& = PSL 

BICW #*C<PSLSM_IV>, R4, -(SP); Save current IV bit on stack 

BICPSW #PSLSM_IV ; Turn off integer overflow trap 

EMODH #1, #0, (RS), R4, -(SP) ; R4& = low 32 integer bits of 'x:/45 
: (SP) = fractional part of ‘xi/45 

SUBH2 (SP)+, (R5) ; (R5) = Integer part of :xi/45 = 1 

BISPSW (SP)+ 3; Restore IV bit 

BLBC R4, a) 3 

ADDH2 4 ; (RS) = 1 +1 

BRB 


NO_OVERFLOW: 
EMODH 4H_1_0V_45, #X_1_0V_45, RO, R4, (R5) 


SSS DPA DADA AAA DM BB BB BEE EAN 


FE RM OODNAUES WIN 8 OD NAUE WIN CO OD NAUSEA SO OONAUE WN 00 


3; R4 = I = integer part of 'X:/45 

7 LBC 4, CvT : Branch if ogtant its are even 
75 ADDL3 #1, R4, (RS) ; (R5) = 1 + ! 
id: CVTLH = (R5), (Ro) : (R5) = 1 + 
77 ODD MULH H_45, (R5) $ one} s Lesden 
A: SUBH2 RO, (RS) : (R5) = ¥ 
ft Bits @*xXFB, R4 ; Save only Last three octant bits 

§ CVT: CVTLH R4, (RS : (RS) = 1 

EVEN: Lng H_M45, (RS) ® nat = -45] 
4 ADDH 0, ik : (R5) = Y 
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54 F8 BF 4 i i pice #*XFB, RG 3; Save only last three octant bits 
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dD 10 
MTHSHSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:39: AX/VMS Macro V04-00 Page 44 
2-007 bEGREE. POLYNOMIALS -SEP=1984 rat 38 £4 MTHRTL. SRCIMTHHSINCO.MAR;1 (26) 
1453 1905 ~SBTTL DEGREE_POLYNOMIALS 
1288 1909 | 
1453 19 8 P_COS_D 
65 EC68 CF 71FD 1453 190 CMPH 4H_90_OV_PI, (R5) ; Compare 90/pi with Y | 
41 18 1459 1910 GEQ 2$ 3 Doub e aPhS5'sion isn't needed 
50 65 65 65FD 1458 1911 MULH3 (R5), (RS), RO : RO/R 
5 DD 146 1316 PUSHL R5 : Save_pointer 
F327 CF gp 0 75FD 1462 191 POLYH RO, #COSDLN1, COSDTB1 ; RO/R3 = Q(Y*2) 
32 88 mre jeer 1918 i nt Sa Bo a adel | 
08 Ad 7¢ 14975 1316 CLRQ_ (RSS : (R5) = YHI 
7E 6€ 65 63FD 1473 191 SUBH3 (R5), (SP), =(SP) : (SP) = YLO 
10 AE 69 60FD 1478 1918 ADDH2 (R5). 16(SP) : 16(SP) = Y + YHI 
&E 5 64rd 147D 1919 MULH2 (SP)+, (SP) : (SP) = YLO#(Y + YHI) = A2 
03 13 1481 1920 BEQL 1$ : Check for A2 = 0 
6f oD A2 1483 1921 SUBW 3 ss#13, (SP) : (SP) = A2/2*13 
5 : 62FD 1486 19 g 1$: SUBH2 (SP5+, RO ? RO/R3 = O(Y*2) = A2/2*13 
65 65 64FD 148A 19 MULH2 (R5), (R5) : (RS) = WHI eg 
65 OD A2 148E 1924 SUBW #13, (R5) + (R5) = YHI*2/2°13 
65 08 62FD 1491 1925 SUBH2 #1, (R5) + (R5) = =(1 = YHI*2/2°13) 
50 65 62FD 1495 1926 SUBH2 (R5), RO + RO/R3 = HCOS(Y) 
QOAD 31 1499 19 f BRW CHECK_DEG 
50 65 65 65FD 149¢ 1998 2$: MULH3 (R5), (R5), RO : RO/R3 = Y*2 
Bo By Ath t Heeler” ° 
Py ave ter 
F214 CF OC 50 75FD = 14A5 193¢ POLYH RO #COSDLN2, cospTB2) ; ROVR3 = aesscy) 
55 8B 9 14AC 193 MOVL (SP)+ $ Restore pointer 
009 1 14Ar 1934 BRW ChECK “oe 
50 08 70FD 1283 1936 3S: MOVH #1, _ RO ; RO/R3 = HCOS(Y) 
0090 31 1486 1937 BRW CHECK_DEG 
1489 1938 
1489 1320 N_COS_D 
65 €CO2 CF 71FD 1489 1941 ~ ~ CMPH H_90_OV_PI, (R5) ; Compare 90/pi with Y 
rt: 18 14BF 13e8 BGEQ 2s ; Double precision isn't needed 
30 $2 OED lees 1904 ee i ioce eatnuns 
3 Vv 
F2C1 CF Ob 50 75FD 1268 1945 POLYH RO, #COSDLN1, COSDTB1 : ROVR3 = aCY"2) 
7 88 rRD 14be 1969 1 A 3 a s Pegpeee pene 
08 A? 7C 1406 1948 CLRQ_ 8 (RSJ : (RS) = YHI 
7E 6€ 65 63FD 1409 194 SUBH3 (R5), (SP), -(SP) : (SP) = YLo 
10 AE 65 60FD 1406 1950 ADDH2 (R5). 16(SP) : 16(SP) = Y + YHI 
6E 3 64rd 14E3 1951 MULH2 (SP)+, (SP) : (SP) = YLO*CY + YHI) = A2 
13° 147 1958 BEQL '1$ : Check for A2 = 0 
g D A2 1469 195 SUBW 0S ssa#'13, (SP) : (SP) = A2/2°13 
t 62FD 14e¢ 1954 1$: SUBH2 (SPS+, RO t RO/R3 = Q(Y*2) = A2/2°13 
65 65 64FD 14F0 1955 MULH (R5), (R5) : (R5) = WHI Tg 
65 OD A2 14F4 1956 SUBW #13, (R5) + (R5) = YHI*2/2°13 
§ 8 FD 14F7 19 SUBH2 #1, (R5) > (R5) = =(1 = YHI42/2413) 
5 FD 14F 1958 SUBH3 RO, (R5), RO + RO/R3 = -HCOS(Y) 
1 130 1959 BRW CHECK_DEG 
135 1961 2$ MULH3 (R5), (RS), RO : RO/R3 = Y*2 


an ZF 


E 10 
MTHSHSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 01:39:11 VAX/VMS Macro v04-00 Page 45 
2-007 BEGREE POLYNOMIALS aot 9453833 MTHRTL.SRCJMTHHSINCO.MAR; 1 ° (26) 
14 13 1508 1366 EQL 3 : Check for Y = 0 
3 DD 150A 19 PUSHL R ; Save_pointer 
F1AD gf C 75FD 1 1964 POLYH RO, PEOSDLNE, COSDTB2 é RO/RS = HCOS(Y) 
0 zoo F af 1513 1965 XORW #*x8000, RO : RO/R3 = -HCOS(Y) 
t p 151 1966 MOVL  (SP)+, RS : Restore pointer 
0028 1 : 18 139 BRW CHECK_DEG 
50 €B3D CF 70FD 151€ 1369 3$ MOVH H_M1, RO 3; RO/R3 = HCOS(Y) 
0022 31 1524 1970 BRW CRECK_DEG 
1 7 1971 
1527 1376 N_SIN_D: 
65 65 72FD 1527 197 MNEGH (R5), (RS) ; (RS) = -¥ 
3 B 1974 P_SIN_D: 
50 65 65 65FD 1528 1975 MULH3 (R5), (R5), RO 3 RO/R3 = Y*2 
oe Ce 15 8 1378 BEQL CHECK_DEG 
os Oe Fa 197 PUSHL R5 ; Save _pointer 
F335 CF OD 0 75FD 1534 1378 POLYH RO, #SINDLN, SINDTB : RO/R3 = PC(Y*2) 
55 aE DO 1538 197 MOVL (SP)+, RS : Restore pointer 
50 65 64FD 153 1389 MULH2 (R5), RO 3 RO/R3 = YeP(Y¥*2) 
65 06 A2 154 1981 SUBW & (R5) $ ) = ¥/2* 
50 65 60FD 154 198¢ ADDH2 (R5), RO : 
1549 198 
1549 1984 
1549 1985 CHECK_DEG: 
5E 55 D1 1549 1986 CMPL R5, SP 3; If RS and SP are equal than the 
03 12 154C 1987 BNEQ 1$ ; =the polynentat was invoked via 
5E 10 CO 154€ 1988 ADDL #16, SP : a branch instruction an the 
05 1551 1989 1$: RSB r stack should be cleared 
1386 1990 
1552 1991 


(R5) = Y 
: RO/R3 = HSIN(Y) 


50. (08 7ore 


50 €B804 CF 70FD 
05 


52 OC 

00000000 ' GF 2 Fe 
04 50 4 

52 04 AD C 
5004 

0D 52 06 €1 

E 0D 

7E O0°8F 9A 
Q0000000'GF 02 fe 


o 
— 
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=< 
zv 
> ted 


GN 
or 
mm 
qo" 
ze 
mo 


: Underflow; if 


MOVPSL Re 
cits #0, G*MTHSSJACKET_TST 


BLB RO, 10$ 
MOVZWL SFS$W_SAVE_PSW(FP), R2 
CLRL = RO 


— ed ee a ts a SS mao 
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COCD 08 NINN NSN NP AAA HII 
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AIF WN $$ O DOONAN EW OOD “IOUS UW OOONOUSW 3 
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-SBTTL DEGENERATE SOLUTIONS 


MOVH #1, RO 
RSB 


MOVH  4H_M1, RO 
RSB 


BBC #6, R2, 208 
PUSHL (SP) 
MOVZBL #MTHSK_FLOUNDMAT, -(SP) 


CALLS #2, G*MTHSSSIGNAL 
RSB 


- END 


user has FU set, signal 


1:39:11 AX/VMS M 
1:38:20 UMTMRTL. Ske 


3; Answer is 1 


; Answer is -1 


ro V04-00 
JMTHHSINCO.MAR; 1 


error. Always return 0.0 


: 2 = user’ 


urn 


R s 
RO = TRUE if 
branch if user di 
et user PSL saved by CALL 

0 = resuit. LIBSSIGNAL will save in 
air anata so any handler can 


xup 
; has user enabled yeaee tne under flow? 
yes, return PC from specia 

trap code for hardware float 
; underflow convert to MTHS$_ 
(32-bit VAX-11 
: pt a (condition, PC 

re 


acket routine'’s PSL 


ar peace routine 


l routine 


Page 4 


- 
F LOUNDMAT 
exception code) 
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Macro Page 48 


H 10 
MTHSHSINCOS ; Floating Point Sine, Cosine and Sincos a eb 9 01:39:4) AX/VMS ro V04-00 
-SEP-1984 11:25:20 CMTHRTL.SRCIMTHHSINCO.MAR; 1 (28) 


Symbol table 


SINDTB 4 fh 02 
SINE = 
SINLENC 
SINLENR = 
SINTBC § 2 
SINTBR R 
LL_COS C08 R 
SMALL_COSD D&C R 
SMALL “SI g R 0 
SMALL_SINCOS ASS R 0 
SMALL_SINCOSD coe R 
SMALL_SIND 8 D $ R 
SUBTRACT op R 
UNFL 3 1 R 0 
W_ADJUST = 9 3FF 
W_MAX_ WEIGHT = 00 $08 
W_TERA_ WEIGHT = 000000 
W-YLO_GEIGHT = 00000040 
x = 44 40 
X_1_0V_45 = 0000C16C 
poem em emer ema woa = + 
! Psect synopsis ! 
¢wmoeocan aan maonranoeca + 
PSECT name Allocation PSECT No. Attributes 
» OBS « 00000000 ( 0.) 00 ¢ O.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
SABSS 00000000 ( 0.) O1¢ 14.) NOPIC USR CON ABS LCL NOSHR EXE RD WRT NOVEC BYTE 
_MTHSCODE 00001582 ( 5506.) 02 ¢ 2.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
$e omeeceee sec eee cen eae ese 
; Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 32 0:00:00.11 0:00:01.05 
Command processing 114 B08: 80 8) B28: 3.68 
Pass 235 0:00:06.90 0: 8: 8-88 
Symbol table sort 0 00:00: 2-39 0:00:00. 
Pass 2 400 0:00:04. 0:00:14.6 
Symbol table output 19 :00:00.1 0:00:00.1 
Psect synopsis output :00: 8 3: :00.04 
Cross-reference output :00:00. :00: 08 
Assembler run totals 80 00:00:12.6 0:00:40. 


The working set Limit was 1050 pages. 

46029 bytes (90 pages) of virtual memory were used to buffer the intermediate code. 

There were 20 pages of symbol table space allocated to hold 200 non-local and 68 local symbols. 
source Lines were read in Pass 1, producing 40 object records in Pass 2. 

10 pages of virtual memory were used to define 9 macros. 
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1.10 
MTHSHSINCOS : Floating Point Sine, Cosine and Sincos 16-SEP- 1: AX/VMS v04-00 Page 49 ef 
VAX=-11 Macro Run Statistics ° g-8 ety 9 +33: 4 YRTARTE eRe STH HHSINCO.MAR; 1 (28) 


$emweeeeeeeceeoeccccoccccce + 
! ; Macro Library statistics ; 

Macro Library name Macros defined 

-$255$DUA28 : CSYSLIBISTARLET.MLB; 2 5 

131 GETS were required to define 5 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHHSINCO/OBJ=OBJ$:MTHHSINCO MSRC$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) +MS 
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